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INTRODUCTION 


A SINGLE crystal of CsI (TI) supplied by Harshaw Chemical Co., was studied 
at room temperature for its luminescence spectra under ultra-violet excita- 
tion and scintillation spectra under y-rays. As pointed out by Sciver W. J. 
Van,° the identity of luminescence centres in the two types of excitation 
cannot be assumed unless it is established that the emission spectrum in 
fluorescence under ultra-violet and that in scintillation under y-rays are the 
same. There does not seem to be any systematic spectroscopic study of 
the scintillations excited by y-rays in Thallium-activated alkali halide crystals, 
though these crystals are now in constant use for y-ray scintillation spectro- 
metry. The purpose of these studies was to compare the emission spectra 
of CsI (Tl) crystal under ultra-violet and y-ray excitations. If the emission 
spectrum of scintillations under y-rays had more than one band, it would 
be of great interest to determine the decay time and decay characteristics 
of each of the emission peaks separately. For this reason, it is necessary 
to measure the emission spectrum of scintillations in some detail. For 
purposes of comparison, the emission spectra of Thallium-activated CsBr, 
Nal and KI were also studied. 


EXPERIMENTAL ARRANGEMENT 


A Steinheil 3-prism spectrograph with collimator and camera having 
apertures ~f/3 and ~/f/4 respectively was employed and the spectra 
in the visible range 4000 A-6000 A were photographed on a Kodak Tri-X 
fim. The average reciprocal dispersion of the spectrograph in the visible 
range was ~50 A/mm. A condenser with an aperture ~f/1-5 collected 
light from the crystal on to the slit of the spectrograph. As the spectra 
at room temperature were expected to be broad continua and the light 
output was very low especially with y-excitation, the slit-width was kept as 
large as 4001. Co ( ~7 m.c.) provided a y-ray source while suitable filter 
solutions were employed to cut off visible radiation and transmit mainly 
ultra-violet from an intense high-pressure quartz mercury lamp. The crystal 
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was kept in a light-tight housing as exposures of the order of a few days were 
necessary to obtain spectra under y-ray excitatior. 


The liquid filters used with the, ultra-violet source were as follows: 


(1) (CoSO, + NiSO,) solution in water in a quartz cell having a path- 
length of approximately 2cm. transmits the radiation between 3600A and 
2300 A with a peak transmission of ~ 75% at 3000A. This filter could 
thus entirely remove the visible continuum in the spectrum of Hg-source 
and reduce considerably the intensities of lines in the visible region. This 
was necessary to prevent visible light reflected and scattered from the 
crystal, masking the genuine luminescence spectrum. 


(2) (I, + CCl,) in a similar quartz cell transmits between 4400A 
and 3000A with a peak transmission of ~ 60% at 3800 A. 


These filters were employed one at a time between the Hg-source and 
the crystal. With them, the space in the visible range was made clear, as 
it were, for the true luminescence spectrum to appear. Duration of 
exposures required were of the order of a few minutes. 


RESULTS 


The spectrograms thus obtained were microphotometered on a Hilger 
non-recording instrument. The results are given in Fig. 1(a)-1 (g). 


The scintillation spectrum of CsI (Tl) excited by y-rays is continuous 
between 4000 A and 6000A with at least three peaks at approximately 
5950-50 A, 5700-50 A and 4750-50 A. The three peaks are very well 
pronounced. There is probably a fourth peak at 4400 A but this is rather 
doubtful since it appears as a small shoulder to 4750 A band [Fig. 1 (a)]. 
The luminescence spectrum of CsI (TI) is described as ‘white’ by J. B. Birks® 
while the range is given by Hahn and Rossel* as 4200-5700 A. 


When fluorescence is excited by ultra-violet light with a continuous 
range between 2300 A and 3600 A, the emission spectrum is continuous with 
peaks at identical wavelengths as in scintillation spectrum. These peaks, 
however, are not as well resolved as in y-ray excitation experiment [Fig. 
1(b)]. With the exciting spectrum comprising of wavelengths between 
3000 A and 4400 A, only a part of the emission spectrum appears, i.e., only 
the wavelengths greater than 5600A are excited. Nevertheless, in this 
region the yellow bands appear again at identical positions as with y-ray 
or deep ultra-violet excitation [Fig. | (c)]. These two peaks are certainly 
more prominent than in Fig. | (0). 
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CONCLUSIONS 


The fluorescence and scintillation spectra for CsI (TI) show emission 
bands centred around the same wavelengths. This experimental observa- 
tion indicates that the fluorescence and scintillation centres are one and 
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the same. The same fact was also observed for Nal(TI) and KI(TI) 
crystals. 


Comparison of Fig. | (b) and | (c) shows that only a part, i.e., the longer 
wavelength portion, of the emission spectrum is excited by ultra-violet 
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light transmitted by I, + CCl,) filter, while the complete spectrum appears 
with (NiSO, + CoSO,). Hence, the luminescence in the range of wave- 
lengths less than 5600 A is excited by the wavelengths between 2300 A and 
3000 A. The experiments with monochromatic ultra-violet excitation would 
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probably yield the wavelength at which 44750A band begins to appear, 
Thus, the transition energy of the excited state responsible for 4750 A band 
‘would be obtained. Some calculations of the shapes of the potential energy 
curves for the ground and the excited states to account for the transition 
energy, so obtained, could then be attempted. 
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As only the yellow bands can be excited by wavelength greater than 
3000 A while the 44750 A band appears when the shorter wavelengths are 
present in the exciting spectrum, it can be suggested that the excited state 
for the yellow emission lies in between that for 4750 A emission and the 
ground state. It is, of course, assumed that all the optical transitions belong 
only to one type of luminescence centre, presumably Tl* substituting Cs*. 
The transition energies corresponding to 5800 A and 4700 A emission are 
roughly 2-1 ev. and 2-6ev. respectively. The absorption data for CsI (TI) 
as given by Pringsheim! consists of three peaks centred at 2990 A, 2690 A 
and 2410 A. The transition energies corresponding to these three peaks 
are approximately 4-14 ev., 4-6ev. and 5-14 ev. respectively. Our experi- 
ments indicate that the absorption of 2990 A is responsible for the lumine- 
scence emission in the yellow region, with about 2ev. dissipated in the 
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crystal corresponding to one yellow photon emission. The correlation bet- 
ween 14700 A emission and the absorption can be established only with 
the monochromatic ultra-violet excitation experiments. 


When Fig. | (a) is compared with Fig. | (d), we find that the emission 
spectrum of CsI (TI) shows a remarkable similarity to that of CsBr (TI), 
although the latter shows a peak at 4500 A instead at 4750 A. Also from 
Fig. | (f) and Fig. 1 (g), we observe that the emission of NaI (TI) is almost 
identical with that of KI(TI) under ultra-violet excitation both showing a 
single band centred at 4300A. Incidentally, this value agrees fairly well 
with those given by Birks as 4100-+-400 A for Nal (TI) and 4100+450 A 
for KI(Tl) and also with the values reported by Milton and Hofstadter 
as 4100-+850 A for NaI (Tl) and 4000-+600 A for KI (TI) under y-ray excita- 
tion from Ra source. 


We note that CsI (Tl) and CsBr(T!) have a b.c.c. lattice while Nal (TI) 
and KI(FI) are f.c.c. With high concentration of Tl, KCl (TI) which is 
also f.c.c. emits a band at 4700 A.2 In spite of these differences in details, 
it appears that the essential structure of the luminescence spectrum depends 
upon the type of crystal lattice or more precisely on the number and dis- 
position of the nearest neighbours of the substitutional TI* ion. The small 
differences in the position of the peaks may be attributed to the types, sym- 
metry and distances of perturbing ions. 


SUMMARY 


Fluorescence and scintillation spectra of a single crystal of CsI (TI) 
at room temperature excited by ultra-violet (with A = 2300 A-3600 A) and 
y-rays from Co respectively are presented. They are continuous in the 
visible region (4000 A-6000 A) with at least three well-pronounced peaks at 
5950-450 A, 5700-+-50A and 4750-+-50 A. There is an indication of a fourth 
peak at 4400 A. The same peaks appear in both types of excitation. Hence 
the luminescence centres in the two processes seem to be identical. With 
the exciting ultra-violet comprising of A> 3000 A, only the two peaks in 
yellow region appear. Hence, it seems that the excited level responsible 
for the yellow emission lies in between the ground and the excited state 
responsible for 4750 A emission. Luminescence spectra for CsBr (TI), 
Nal (Tl) and KI(TI) are also presented. It is concluded from the com- 
parison of all the spectra obtained that the type of emission spectrum depends 
on the crystal structure. 
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ABSTRACT 


The pre-sunrise phenomenon of a rapid drop in the maximum ionic 
density of F-layer followed by an increase in the semi-thickness and 
height of the layer and stratification is described. The phenomenon, 
which is found to occur prominently during periods of low or moderate 
solar activity and during winter months, is discussed in relation to the 
distortion of the overhead Sq. current system and vertical drift of 
ionization. 


1. INTRODUCTION 


VERTICAL soundings of the ionosphere using a N.B.S. type C-3 Ionosonde 
were started at Kodaikanal (Latitude, 10° 14’ N., Longitude 77° 28’ E.) early 
in 1952. The early morning phenomenon of F-layer stratification was 
observed soon after and was reported, on the basis of preliminary observa- 
tions, by one of the present authors (1952). Similar observations have 
been made by Bandyopadhyaya (1959) at Haringhata near Calcutta and 
have been discussed for possible bearing on the E-layer cusps observed there. 
Round-the-clock observations at Kodaikanal were commenced in September 
1955 and intensive observations were made during the 3-hour period pre- 
ceding sunrise in order to determine the nature of the stratification which 
was frequently observed during periods of low solar activity and during 
winter. 


2. EFFECT OF VERTICAL DRIFT ON ELECTRON DENSITY 


The basic continuity equation of electron density N in a simple solar- 
controlled region is 


= = q(t) — aN? — div. (Nv) (1) 
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where q(t) represents the ion production term and aNP is the electron dis- 
sipation rate, p being the exponent to indicate whether the dissipation pro- 
cess is one of recombination (p = 2) or of attachment (p = 1). The term 
div. (Nv) represents the electron transport rate. 


If we restrict our attention to the peak of the layer, Equation (1) reduces 
to 


IN dw 
57 = 4 (1) — aN? — Nm 5 - (2) 


Martyn (1947 a, b, 1948) had indicated that vertical movements may cause 
important perturbations in an ionized region. Appleton and Lyon (1955) 
have discussed the effects of vertical drift in distorting the electron density 
profile of a Chapman-like region. They have shown that the fractional 
change in the level of maximum electron density 54 (Nm) and in its magni- 
tude SN», due to vertical drift are given by 


5h (Nm) = 7a Nz (recombination) 


5h (Nm) = 2 (attachment) 


dNm = l ; ov 1 v 2 ; 
Nm 2aNmn 2 4H? (sn) (recombination) 
5N l ov l v\2 

x. at B ‘Th 2HeE ( 4) (attachment) 


where a and § are recombination and attachment coefficients respectively, 
v is the vertical drift velocity reckoned positive upwards, (dv/dh) is its vertical 
gradient and H is the scale height of the ionized constituent. Equations 
(3) and (4) indicate that the influence of vertical drift on Nm and its height 
is determined by drift characteristics and by the quantity (1/2aN) or (1/f). 


(4) 


3. CHANGES IN F-LAYER BEFORE SUNRISE 


As stated previously, the early morning hours of winter months, parti- 
cularly during periods of low solar activity, are characterised by rapid changes 
of maximum ionic density, height and layer shape. A typical sequence of 
ionograms observed between 0430 and 0630 hours on the morning of 
January 4, 1960, is given in Fig. | (Plate VII). The layer is normal ‘till about 0430 
hours. Shortly afterwards, the ionic density begins to drop and the virtual 
height registers an increase. Around 0600 hours the critical frequency of 
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Virtual Height (Km) 


15 20 25 30 35 40 45 50 55 60 
Frequency (mc./Seéc) 


Fic. 2. Ordinary component of the h’-f traces between 0400 hr. and 0630hr. IL.S.T. on 
January 4, 1960. 


about 2-5 Mc./sec. is registered and rapid changes in the structure of the 
layer are noticed. The layer is stratified and the two strata become well 
defined. Thereafter, the maximum electron densities in both the layers 
increase and the two layers merge. The layer becomes normal by about 
0630 hours. These changes are shown in Fig. 2 where the ordinary com- 
ponents of the h’-f traces between 0400 and 0630 hours are reproduced. 


4. THE DISTRIBUTION OF ELECTRON DENSITY AS A FUNCTION OF 
TRUE HEIGHT 


In order to examine the variations of electron densities at several levels, 
the h’-f records during the interval 0330 hours to 0630 hours for two typica| 
days, January 12, 1956 and January 4, 1960, were reduced to (N, /); profiles 
using Tables of Coefficients given by Schmerling and Ventrice (1959). The 
profiles for January 12, 1956 and January 4, 1960, are given in Fig. 3 a and 
Fig. 3 b respectively. From the (N, /)¢ profiles for one of these days, viz., 
January 4, 1960, the height variations of constant ionization levels were 
deduced and the resulting curves are presented in Fig. 4. It will be seen 
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Fic. 3(a). (N, A), profiles between 0300 hr, and 0700 hr. I.S.T., deduced from h’-f records 
of January 12, 1956. 


from the curves of Fig. 3 a and Fig. 36 that considerable variations in 
height and thickness occur near the layer peak from 0430 hours and these 
variations are maximum around 0600 hours. Fig. 4 shows some interest- 
ing features. Perturbations in the heights at different ionization levels take 
place. At 0600 hours when the layer gets stratified, the variations in height 
appear to be maximum. It is also noticed that the layer appears to rise 
as a whole with a velocity of about 100 km./hr. Simultaneously, the critical 
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Fic. 3(b). (N, /); profiles between 0330 hr. and 0630 hr. L.S.T., deduced from h’-f records 
of January 4, 1960. 
frequency drops to a low value of about 2-5 Mc./sec. at the time of strati- 
fication. 
5. DISCUSSION 


(i) Changes in the layer shape due to vertical movements.—As mentioned 
earlier, the effect of vertical drift velocity is to perturb the ionized layer, the 
electron density changes being represented by the continuity equation (1). 
It will be further seen from the Appleton and Lyon equation (4) that irres- 
pective of its sign, the effect of drift velocity v will be to spread out the ioniza- 
tion, The rise or fall in height of peak ionic density will depend upon the 
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Fic. 4. Height variations of constant ionization levels of January 4, 1960, during pre-sunrise 
F-region “‘splitting’’. 


sign of v, equation (3). The magnitude and sign of the gradient of drift 
velocity will determine the overall effect of distortion on peak ionic density. 


Using the electron density-height profiles the time variation of the 
heights of peak ionic density was obtained. It was found that the peak 
rises rapidly from 270 km. at about 0430 hours to over 450 km. at 0600 hours. 


According to Martyn (1955), the effect of height gradient of electron 
loss (da/dz) is to move the level of Zm upwards. It can be seen from 
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the constant ionization levels (Fig. 4) that just before sunrise, the level of 
Zm is very much lifted up, indicating thereby, that the magnitude of (da/dz) 
is very much increased. 


Following Martyn, we shall write 
a = ae (5) 
Differentiating the above, we get 


da a 
hl Se (6) 


If the magnitude of (da/dz) is large, the value of “a” should be small. This 
can be examined in a different way. Following Marasigan (1958), we may 
write that at the point of inflexion 


© (aNP) = 0. (7) 


Differentiating (7) we get 


da a 92 
az ~~ NP * 9z‘NY). 


From (6) and (8) we get 
N 


(32) 


a= 


= 1, attachment process). 


=a (p 
CG: 


It is an observed fact that just before sunrise the magnitude of N is very 
much reduced while (9N/dz) is much increased, thereby indicating a small 
value for “‘a’’ just before splitting. It appears, therefore, that the large upward 
shift in the level of maximum ionization is followed by splitting of the layer 
due to large value of height gradient of loss coefficient (da/2z). 


(ii) Magnitudes of vertical drift velocity v and its gradient (dv/dz).— 
Using the relation 54(Nm) = (v/8) where £8 is the attachment coefficient, 
values of » were computed at 10 minutes intervals on January 4, 1960. 
Substituting these values of v in the continuity equation, the magnitude 
and sign of (dv/dz) at the height of maximum ionization at the same intervals 
were also computed. These values are given in Tables I a and Ib, 
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TABLE I (a) 





Magnitude of drift velocity 
gradient (dv/dz) on 
January 4, 1960, in sec.-! 
at 10-minute interval 





105 x 
0400 —14°9 
10 —13-°5 
20 —11+5 
30 —13-1 
40 —19+4 
50 —18+4 
0500 —16°8 
10 —15+2 
20 —13°5 
30 —11°7 
40 — 9°5 
50 — 61 
0600 — 2-9 





It will be seen that starting from about 0440 hours v was positive and 
its magnitude increased rapidly. At the same time (2v/dz) approached zero, 
starting from a negative value, thereby indicating that the second term of 
the right side of the equation 


8Nm _ __ | 2? — an(s) 
Nm B22 2H*\B 


was predominant in determining the electron density changes at the peak. 


(iii) Distortion of the Sq. current system and the associated movements.— 
Martyn and others (1952) have indicated that the morphology of the F- 
region may be profoundly affected by the events in the E-region becaus¢ 
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TABLE I (5) 

Magnitude of drift velocity 
Time v in units of 10? m. 
LS.T. sec.-! on January 4, 1960, 


at 10-minute interval 





0400 —1-95 
10 —1-45 
20 —0-80 
30 +0-60 
40 +0-95 
50 41-27 

0500 +1-42 
10 +1-50 
20 +1-66 
30 +2-00 
40 +2-60 
50 +4-40 

0600 43-52 





the electric fields would be communicated to F-region along the lines of 
magnetic force. This is supported by the experiments of Maple, Bowen 
and Singer (1952) and more recently by those of Van Allen and Cahill (1958) 
which indicate that the E-region is the seat of strong electric currents and 
that these currents are responsible for the diurnal quiet-day Sq. variations 
of the magnetic elements. Any change in the horizontal intensity at low 
latitude stations, therefore, indicates a distortion of the Sq. current inten- 
sity, the magnetic field being related to the electric field by 


4H = 2/11 = 211K4E. (9) 


Further, the magnitude of the drift velocity v is determined by the east-west 
component of the Sq. current and the strength of the horizontal component 
of the earth’s field, 
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In order to examine the relationship between the abnormal rise in the 
peak of electron density in the pre-dawn period and any distortion in the 
Sq. current intensity, the following analysis was undertaken. Thirty days 
in January, 1960, for which complete data were available were grouped into 
two categories depending upon whether the F-layer critical frequency 
immediately preceding sunrise was below or above an arbitrarily chosen 
value of 3-0 Mc./sec. There were nine days of the first category and twenty- 
one days of the second. The departures 4foF2 from mean values for these 
days were derived at half hour intervals between 0300 and 0700 hours LS.T. 
Likewise, the departure 4H in gammas in horizontal force recorded at 
Kodaikanal were also obtained for the two groups of days at half hour 
intervals. These departures 4foF, and 4H are plotted in Fig. 5. It can 
be seen at once that on nine days when foF, is negative, 4H is also negative, 
the maximum dip in 4foF2 occurring at about 0600 hours, about an hour 
and a half later than the dip in 4H. It is interesting to note that it is just 
about 0430 hours when the maximum dip occurs in H, that the height-changes 
in F-layer begin to take place (Figs. 4 and 5) presumably due to distortion 
of the Sq. current intensity in the dynamo region. The maximum departure 
in foF2 occurs after an interval of about an hour and a half which is very 
close to the relaxation time [7 = (1/2aNm)] where a is the effective loss 
coefficient taking into account the vertical drift 


1 ww 
(. =%T yh) * 


It, therefore, appears that the abnormal rise of the layer in the pre-dawn 
period is associated with the additional movement in the region caused by 
the distortion in the Sq. current intensity. 


(iv) Influence of solar activity on the occurrence of the pre-sunrise pheno- 
menon.—The pre-dawn increase in the height of the F-layer and the associated 
splitting were rarely seen during 1957-59, a period of high solar activity. 
It is known that the rate of ion production is in general proportional to the 
intensity of solar ionizing radiation. The magnitude of the first two terms 
on the right side of the continuity equation (1) will be large as compared to 
the last term and the effect of the transport term div. (Nv) is likely to be 
small. Further, at sunspot maximum there may exist (Allen, 1953) a high 
resistance to movements due to electric or magnetic damping by the higher 
electric conductivity. The contribution of the movements to the rate of 
change of ionization becomes progressively smaller with increasing solar 
activity while the residual overnight ionization densities are comparatively 
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Fic. 5. Upper curves: Variation of A foF2 for January 1960: (a) for 21 days when pre- 
sunrise minimum of foF2 was >3-0 Mc./sec. and (5) for 9 days when the parameter was fcund 
to be < 3-0 Mc./sec. 

Lower curves: Variation of AHat Kodaikanal for January 1960, averaged for the same 
days as in the upper curves (a) and (6). 


much larger. The conditions favourable for bifurcation are, therefore, 
not fulfilled during the periods of high solar activity. 
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SoME years ago Otto Szdsz! showed that the Bessel function J, (x) satisfies 
the inequality 


An (x) = Un GO — Ins @) Ina) > op Un OP, 


n>0, x real (1) 
wherefrom he deduced the interesting inequality 
An€x)>0 for n>0O, x real. (2) 


The Turan inequality (2) for the function Jn (x) is an immediate consequence 
of the identity? 


co 


An (x) = 4 pM (n + 2k + 1) (Sniek+s (x))? (3) 


k=0 


which represents A» (x) as a series of positive quantities whenever n > — I. 
V.R. Thiruvenkatachar and T. S. Nanjundiah* obtained a new positive repre- 
sentation for An(x) by proving the identity 


An (2) = 17 (In)? + 4-3 Unie )? 


+203 [n+k—-YO@+k+ DP Crue @)? 


k=2 


and the above result (1) of Otto Szasz follows from this immediately. 


While the above representations for An (x) are quite interesting we may 
also proceed to deduce the inequality (2) by considering the behaviour of 
the function 4, (x) itself. This procedure leads to several questions con- 
cerning the sequences of relative maxima and minima of A» (x) and the 


* Present address: Regional Engineering College, Warangal (A.P.), India. 
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aim of the present note is the study of these sequences. We prove in parti- 
cular (Theorem 4 below) that the r-th relative maximum of A» (x) is larger 
than the r-th relative maximum of /Ayn,, (x), for any fixed value of r. This 
answers the question first raised by John Todd for the classical orthogonal 
polynomials and settled for the Hermite function Hy (x) by Otto Szdsz* and 
for the Turan expression for the Hermite function by the author. 


THEOREM |: When n> 0, the relative maxima of An (x) occur at 
the zeros of Jn_, (x) and the relative minima occur at the zeros of Jn,, (x). 


On differentiating the expression 
An (x) = (In (x))? — Ins (X) Ina () (5) 


twice, we get 


£ An (x) = : Inia (x) Ini (x) 


a 2 d 
dx? An (x) = x Inia (x) a In-1 (x) + Stes (x) = Inia (x) 


— 3 Jan @) Ina. 0 


Let 0 < x1,5 < X2,5 < ... denote the zeros of J; (x) in the ascending order. 
From (7) we see that 


(Se An ©), ae ms ae {(; In )) } L=Lr,n-1 lie 


(i An (*) beitaisias, 4n {(; In ©) bene = > 0. 


Hence the theorem. 


The statement in the above theorem is reversed if »<0. In what 
follows we assume that n> 0. 


If we denote the successive relative maxima of An(x) by Min 
M,,n, -.- and the successive relative minima by m,n, mz,n, ... respectively, 
we have 


Mr,n = An (Xr,n-1) = (In (Xr, n-1))?, (8) 








Iti- 


ger 
his 


ind 
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be 
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and 
Mrn = An (Xr,n+1) = (Sn (Xr,ni1))?- (9) 
We may also observe incidentally that 


Mr,n-1 = (Sn (xr,n))? - (Snir (xr,n))? = My, nit. 


Since the quantities M,y,, and my,» are positive for all values of 7, it 
follows that 


An (x) > 0 (10) 


which is Turan’s inequality for the Bessel function J» (x). 


THEOREM 2: The sequences of relative maxima and relative minima 
of An(x) are decreasing beyond a certain value of r. More specifically 


Mryn > Mein if Xrn1 > € = V2n(n — 2) 
and 
Mn > Mian if Xp,n > 7 = V2n(n + 2). 
If 0<n <2, we define = 0. 
Consider the function 
f (8) = An (=) + 22 Ina) 


Then, by differentiation we get 


f' (x)= = Mine (x) Ins (x) — $ (Sn (x))? + co Ins (x) Ms In—1 (x) 


‘a a (2n (n — 2) — x*) (Ina (@))? 


so that f(x) is increasing in 0 < x < € and decreasing in x >. Also at 
the relative maxima of An (x), f(x) = An(x) and hence the assertion that 
Men > Mrisn if Xris,n1 > Xr,n1 > To prove the corresponding 
result for the relative minimum, consider the function 


2 ; 
g(x) = An(x) — = (Sn (x))?. 
Differentiating this we have 


#' (x) = 2.2m +2) — x4) Inna)? 
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so that g(x) is decreasing in x > y. Since g(x) coincides with Ax (x) at all 
the relative minima of A,,(x), it follows that my,, > Mrin if Xra,n4 > 
Xr.nti > 1: 


This completes the proof of the theorem. 


THEOREM 3: The r-th relative maximum of Ay(x) is greater than 
the r-th relative minimum, i.e., Mr,» > Mr,n. 


We have An (Xr,n-1) = Mr,n and An (Xr,n4) = Mr,n and (d/dx) An (x) 
= (2/x) Inia (x) Ins (x). From the interlacing properties of the zeros of 


Bessel functions we know that® ? 


Xra,nt < Xr,n-1 < Xr,n < Xen < Xt, n-1- 
Hence in the interval (x7,n-1, Xr,ni1), Sgn (d/dx) An (x) = (—D)” 71 = - 1, 


so that An (x) is decreasing in the interval and An (Xr,n-1) > An (Xr,nu) 
or My,n > M,r,n which proves the theorem. 


By employing the functions f(x) and g(x) already introduced, we can 
show further that 


8n* ; 
M;,n > (1 + _ ) Mn if Xrna> § 


8n* \-1 : } 
Mr.n > (1 ~ 2. ) Mr,n if Xe,n-+4 > 2. 
“*~ T,N-1 

THEOREM 4: For a fixed value of r, the r-th relative maxima of An (x) 
form a sequence of decreasing functions of n, i.e., Mr.n > Mer,ni. We 
have already noticed that Ap (Xr,n-1) = Mr,n. We may see easily that 
Meynvs = Ansa (Xr,n) = An(X7,n). Also from the relation (d/dx) An (x) 
= (2/x) Inu (x) In (x), it follows that in the interval (%p,n-1, Xr,n), 
Sgn (d/dx) An (x) = — 1, so that An (x) is decreasing in the interval and 
An(%r,n-1) > An(%r,n). Whence the theorem. 


Analogously for the relative minima, we have 


THEOREM 5: For a fixed value of r, the r-th relative minima of An (x) 
form a sequence of decreasing functions of n, i.¢., mp,n > Mp,ni1- 


We have noticed earlier that my,n-; = Mr,nyy. Combining this with 
ihe result of the previous theorem, we see that m;,n-; > mr,n which proves 
the assertion. 








l 
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INTRODUCTION 


ORGANIC crystal scintillation spectrometers are commonly used for f-ray 
spectrometry where moderate resolution in energy is sufficient for the purpose. 
However, its utility in y-ray spectrometry is limited by the main mode of 


interaction being the compton effect which gives rise to a continuous pulse 
height distribution. 


A photoelectron spectrum of y-rays from a suitable converter allows 
one to obtain monoenergetic electron lines from which one can determine 
the energy of the y-ray with an accuracy comparable with that of the con- 
version electron line. However, the intensity of y-rays cannot be measured 
with the same accuracy because of the energy dependent angular distribu- 
tion of photoelectrons with respect to the direction of the incident y-ray. 


[t is the purpose of this paper to describe a new possible method in 
which a combination of the photoelectron spectrum and the electron spectrum 
with the organic crystal scintillation spectrometers is employed to measure 
conversion coefficients and in suitable cases, relative intensities of y-rays. 


Such a combination will be called in the following the photoelectron 
scintillation spectrometer. 


The recent report (Ramaswamy, 1950) of a discrepancy between the 
conversion coefficients, derived by the magnetic spectrometric method 
(0-159) and by the scintillation spectrometric method (0-195) for the 279 Kev. 
transition in T1*°*, aroused the interest in checking the possible existence 
of a systematic error in the scintillation spectrometric method compared 
with the magnetic spectrometric one which has often been suspected (Subba 
Rao, 1960). 
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PRINCIPLE OF THE PHOTOELECTRON SCINTILLATION SPECTROMETER 


The distribution of photoelectrons relative to the angle 4 between the 
direction of propagation of the incident photon and the photoelectron 
released is indicated in Fig. 1(b). This distribution is based on Sauter’s 
relativistic equation* (Heitler, 1955) with the modification that in view of 
the rather low energy photons expected to be met with in experiments, 
photoelectron energy is considered in the present derivation to be (hv — E;)/ 
mac? instead of hv/mgc® as in the original one, hv and E; being the photon 
energy and the binding energy of the electron in the i-th shell. Though 
experimental verification is as yet inadequate, the effect of the azimuthal 
($) distribution of photoelectrons is taken on theoretical grounds (Heitler, 
1955) to be negligible in the following. If a suitable converter material 
can be chosen such that hv > 2E;, the probability of emission of photo- 
electrons into the backward hemisphere is very much reduced [Y in Fig. 
1(b); X corresponds to hy <2E;]. It is to be noted that the multiple 
scattering effects which would further distort the pattern of angular distri- 
bution of photoelectrons in the converter do not affect the intensity 
measurements in the exp:zrimental arrangement described below. With 
the restrictions of absorption of photoelectrons in the converter, one can 
employ fairly thick photoelectron converters in order to obtain higher 
photoelectron yield. Further, such a choice (hv > 2Ej) of the converter 
material also results in the photoelectrons having enough energy so that 
absorption losses are minimised. 

sin? 6 13(1— v1 — gp) —26? sin? 9 











1O= (—peosey *2~ ANF T= Boosey” 
in which 
hv — E, hv — Ey hy — Ex 
B iD = + + 
ae 


From these considerations, an arrangement to study the photoelectron 
spectrum was arrived at [Fig. 1 (a)]. It consists of an anthracene crystal 
\" dia. x}" thick with a well of depth of 4”, covered with a 150 » gm./cm.* 
aluminium reflector. The source is placed over a }” thick lead and 1/16” 
copper collimating rings in a perspex source mount with a well-defined 
position for the source holder. The photoelectron converter is placed in 
the well of the crystal. The source-converter geometry is fixed and is such 


*Note.—Sauter’s formula for angular distribution of photoelectrons is taken finally in the 
form: 
A¢ 
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that when the converter is removed electrons can be detected in the same 
geometry. Source-crystal distance is about 7mm. so that the effect of 
electron scattering in air is not very significant for electrons of energy greater 
than 60 Kev. 


In Fig. | are shown the internal conversion spectrum, the photo- 
electron spectrum and the compton background of the 661 Kev. transition 


a wae near ie 26. oe ee 
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detector geometry. Photoelectron spectrum was taken with a 19-2 mgm./ 
cm.” gold-converter placed in the well of the anthracene crystal. Compton 
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background was determined by replacing the copper collimating ring with 
a copper disc of the same thickness (1/16"). Adopting the photoelectric 
absorption coefficients due to White (1953) the total conversion coefficient 
obtained by comparing the areas under the internal conversion line and the 
photoelectron line is a=0-115+ 0-007. Adopting K/L +™M= 4-6, 
ak = 0:092 + 0-01 is obtained. This is in good agreement with the earlier 
measurements by other methods (Strominger, Hollander and Seaborg, 
1958). 


CONVERSION COEFFICIENT OF THE 279 KEV. TRANSITION IN TL?° 


Hg"? produced by the neutron irradiation in APSARA (INDIA) for 
two weeks of spectroscopically pure Hg,O3 was used as the source. The 
following measurements were made 72 days after irradiation. 


In order to minimise the compton background without any loss in 
electron detection efficiency, a thinner anthracene crystal (4” dia. x }” thick) 
with a well (1/16” deep and 3/8” dia.) was employed. The rest of the arrange- 
ment was as described earlier. The internal conversion lines were separated 
from the gross electron spectrum [Fig. 2 (a)] by first subtracting the compton 
background [Fig. 2(6)] under this and then subtracting the continuous f- 
spectrum. With the energy calibration remaining the same in all these 
observations, the true 8-background was obtained by first obtaining the 
B-spectrum in coincidence with the photopeak of the 279 Kev. y-ray and 
then normalizing this coincidence 8-spectrum at the maximum intensity 
point of the gross electron spectrum from which the Compton background 
had been subtracted. The resulting conversion line is indicated in Fig. 2 (a). 
The (L + M....) part of the internal conversion line tails off at the higher 
energy end. This is attributed to the simultaneous detection of 8-particles 
and the conversion electrons in coincidence with them. This is possible 
since the mean life of the excited state in T1?°* is about 5x 10! secs. and 
the overall detector resolution time is about 2x 10~’ secs. In order to take 
into account the conversion electrons removed from the main peak by this 
process, the area under the whole conversion line including the tailing part 
was employed and the total conversion coefficient (a) was determined. 
However, assuming a gaussian distribution for the K-conversion line, 
K/L+M = 2-4 was obtained. This is in good agreement with the value 
of 2:7 obtained by magnetic spectrometric measurements. This indicates 
the correctness of the above procedure. 


The photoelectron spectrum obtained with a 28 mgm./cm.? lead con- 
verter is shown in Fig. 2(6). In Fig. 1 as well as in Fig. 2, all points 
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include statistical errors. Since the photoelectron and internal conversion 
lines were obtained with the same source and in the same geometry, the 
total conversion coefficient a is deduced from the relation, 


_Ne N & wy _ Ne 
Np N’ «, 


a e—_= - € 
© we Np p> 


where N = N’ (same source), «ge = €e’ (same electron detection efficiency 
of nearly 100%), we = w, (same geometry) and Ne and Np are the conversion 
electron and photoelectron line intensities and ep is the photoelectric con- 
version efficiency for lead. Thus a = 0-205 was obtained as an average 


i 





7 ae 


PHOTOELECTRON SPECTRUM 


COUNTING RATE —=— 














PULSE HEIGHT. ——— 
Fic, 2 





Photoelectron Scintillation Spectrometer 249 


of two runs after correcting for absorption of photoelectrons in the lead 
converter and employing the calculated (White) photoelectric cross-sections. 
With K/L+M = 2:7, a2® = 0-15+ 0-02 is obtained. This error is three 
times the mean deviation in the average value so that any possible syste- 
matic errors in subtraction of backgrounds, normalization of f-spectrum 
or uncertainties in the theoretical photoelectric cross-section are taken into 
account. This value is in agreement with the magnetic spectrometric value 
of 0-159 + 0-004. Thus, the photoelectron scintillation spectrometer 
can be useful in the determination of conversion coefficients by a very simple 
and direct procedure. Where uncertainties due to background can be 
reduced still further, quite accurate results may be expected. 


The conversion coefficient of this transition was redetermined by employ- 
ing a 14”X14” Nal (TI) crystal and a point source placed at a distance of 
9-6cm. from the crystal along the axis of the crystal. The source holder 
was held in position by an insulation tape which was fixed at the two ends 
| yard apart. The crystal was provided with a graded shield consisting of 
lead, cadmium and brass and there was no material within 3 yards in front 
of the crystal. Thus background and scattered radiation were reduced 
considerably. The background below X-ray peak was inferred by suppres- 
sing the X-rays by means of a 100 mgm./cm.? tin-absorber placed midway 
between the source and detector [Fig. 3(a)]. This background, consisting 
mainly of compton distribution and the back-scattered peak, is shown by 
dashed lines. Background under the spectrum without the source has been 
subtracted from the spectra [Fig. 3(6)] I & If. After corrections for absorp- 
tion in the material surrounding the crystal and taking into account the 
efficiency of detection (Wolicki and others) a* = 0-153 + 0-012 was obtained. 
Since in this method a* is deduced from the relation 


x — Nx & &x I 


a Se 


> 
y by Qy OK 


with the fairly well-known values of the K-fluorescence yield w, and detec- 
tion efficiencies ex and «, and with the solid angle of detection being the 
same (Q, = Q,) the large deviations could arise mainly through the X- and 
ytay intensities. This can happen if either the background under the 
X-ray peak is uncertain or there is some extraneous contribution to the 
X-ray region. It is not obvious how the same method could lead to results 
different from those of earlier workers by the scintillation spectrometric 
method (Ramaswamy, 1960). This has been repeated three times and all 
of them agree with the value 0-153 + 0-012. 
AS 
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However, it appears relevant to point out that the spectra studied only 







































































a fortnight after irradiation of the source-sample lead to a result aX = 
0-25 + 0-02 in spite of the fact that the Hg!’ activity which contributes 
mostly to this X-ray region had undergone a decay for nearly five half-lives, ( 
This result was obtained by the same method as described above. | 
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From these measurements, it is concluded that scintillation spectro- 
metric methods do not indicate any systematic errors compared with magnetic 
spectrometric method (Nordling etal., Wapstra and Nijgh, 1957) in the 
case of the 279 Kev transition in TI, 7 
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SUMMARY 


An experimental arrangement combining the photoelectron and organic 
crystal scintillation spectrometer methods, called a photoe!ectron-scintil- 
lation spectrometer, is described and its utility in making direct measurements 
of conversion coefficients without a specific knowledge of source-strength, 
source-detector geometry, efficiency of detection and the energy dependent 
angular distribution effects and multiple scattering effects, is described. 


By this method and the conventional method of measuring X- and 
y-ray intensities with Nal (Tl) crystal, the conversion coefficients of the 
279 Kev. transition in T]#°* are found to be a* = 0-15+ 0-02 and a* = 0-153 
+0:012, respectively. These values indicate that there are no systematic 
errors in the scintillation spectrometric methods as compared to the magnetic 
spectrometric results in this case. 


The author is much indebted to Dr. B. V. Thosar for some helpful dis- 
cussions. 
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ABSTRACT 


This paper gives an account of the results of the tests that have 
been carried out to assess the relative importance of the phase angles 
and structure amplitudes in a Fourier synthesis. For this, the phase 
factors, e“*, have been combined in different ways with structure ampli- 
tudes such as when the latter are (a) replaced by the r.m.s. of the structure 
amplitudes (i.e., by 2 f;?), (6) randomly chosen from the set of correct 
structure amplitudes, (c) made us much nti-correlated (about the 
mean) with the correct ones as possible and (d) taken to be those of an 
entirely different structure. The above syntheses which are only modifi- 
cations of the phase synthesis, e“*, have been tested on both the centro- 
symmetric and non-centrosymmetric projections of a known structure. 
They strongly establish the greater importance of the phase angles in a 
Fourier synthesis. 


1. INTRODUCTION 


IN a series of papers from this laboratory (Ramachandran and Raman, 1959; 
Raman, 1959; Raman, 1960; Srinivasan, 1960a; Srinivasan and Aravin- 
dakshan, 1960; hereafter referred to as Parts I, LI, III, [V and V respectively) 
the problem of deconvoluting the Patterson function when a part of the 
structure is known, was considered in detail. In the course of these studies 
it became necessary to analyse the significance of the standard syntheses 
which employ for their Fourier coefficients various standard functions of 
F, such as | F |, e’*, 1/F, etc. It was pointed out in Part I and later proved 
in a more rigorous way in Parts If and IV that the phase synthesis e* 
resembles closely the Fourier synthesis while the modulus synthesis, | F |, 
resembles the Patterson. Though this indicated in a general way the impor- 
tance of the phase angles, a closer study of the relative importance of the 
phase angles and the structure amplitudes was thought worthwhile. The 
results of such a study are presented in this paper. 
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2. DESCRIPTION OF THE VARIOUS SYNTHESES 


The calculations were all made on the centrosymmetric and non-centro- 
symmetric projections of L-tyrosine hydrochloride whose structure has 
been thoroughly analysed in this laboratory (Srinivasan, 1959). The crystal 
data of this compound are: a= 11-1, b=9-1, c= 5:1A, £ = 92°, two 
molecules of C,H,,NO;.HCI per cell, space group P2,. In all the syntheses 
to be described below the phase angles referred to (or the signs in the centro- 
symmetric case) correspond to those calculated for the above structure 
(Srinivasan, 1959) unless otherwise stated. 


(a) The Normalised Phase Synthesis —To start with, a modified form 
of the phase synthesis, e**, was devised. The modification consisted in not 
using unit amplitudes for all reflections as in e**, but a value equal to Vz fj? 
which is the mean for all reflections having the same value of (sin @/A). 
Thus, in effect, the magnitudes of the geometrical structure factors of all 
reflections were made equal, but all of them were multipliea by a mean 
atomic scattering factor. We shall call this synthesis the “normalised phase 
synthesis”, the normalisation being done such that 2 | F |? is the same for 
this synthesis and the Fourier synthesis. The syntheses for the non-centro- 
symmetric and the centrosymmetric projections are shown in Figs. 1 (b) 
and 2(4) respectively along with the actual Fourier syntheses in Figs. 1 (a) 
and 2 (a). 


In carrying out the calculations the value of / fj? for each reflection 
j 


was obtained from a graph. The temperature factor was also applied to 
the amplitudes and it was the same as that employed for the actual crystal. 
It will be noticed that the agreement between the normalised syntheses 
[Figs. 1 (6) and 2 (6)] and the actual Fourier syntheses [Figs. 1 (a) and 2 (a)] 
is extremely good. 


(b) The Random Phase Synthesis ——The interesting question now arises 
as to what would happen if random values of | F | are used along with the 
correct phase angles, but not equal values (apart from the ‘f’ factor variation) 
as in the above synthesis. Such a synthesis which may be called the “random 
synthesis’ was calculated. However, in order to ensure the same statistical 
distribution of structure amplitudes as in the correct structure, the reflections 
were grouped into various ranges of (sin @/A) and in each range the values 
of the actually observed | Fo |’s (Srinivasan, 1959) were randomly permuted 
amongst themselves and these were then used with the correct phase angles 
to calculate the synthesis. The syntheses for the non-centrosymmetric 
and the centrosymmetric cases are given in Figs, 1 (c) and 2 (c) respectively, 
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Surprisingly, even in these random syntheses, the structure is revealed in 
fairly good detail. 
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Fic. 1*. (a) The Fourier synthesis, (b) the normalised phase synthesis, (c) the random 
phase synthesis and (d) the anti-correlated phase synthesis of L-tyrosine hydrochloride projected 
down the c-axis. Contours are at intervals of 2e/A? at chlorine and 1e/A® elsewhere, starting 
from 2e/A?. 


* The different kinds of atoms in the molecule of tyiosine are not marked in all the figures 
but only in Figs. 3 (6) and 4(6) where small dots correspond to carbons big dots to oxygens and 
open circle to nitrogen. 


(c) The Anti-correlated Phase Synthesis —Now, it was thought worth- 
while to go a step further and take the strong reflections weak and vice versa. 
This would be equivalent to making the new set of amplitudes and the 
correct set as much anti-correlated about their mean as possible.* The 
method adopted for getting the maximum anti-correlation was as follows: 
In each range of (sin 9/A) the various reflections were arianged in decreasing 
order of magnitude. The whole sequence was now inverted so that the 


* For a proper explanation of the various statistical terms used in the discussion se¢ 
Appendix I. 
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strongest reflection would become the weakest and so on. However, it 
will not be possible to realise perfect anti-correlation about the mean by 
this process. The possibility of achieving this depends on the nature of 
the probability function and it can be shown that, unless the function is 
symmetrical about the mean, perfect anti-correlation is not possible. In 
our present case, the probability functions for the centrosymmetric and non- 
centrosymmetric cases are not symmetrical about the mean (Ramachandran 
and Srinivasan, 1960) and hence perfect anti-correlation cannot be realised 
for them. However, the maximum possible anti-correlation for any type 
of distribution can be worked out and they are given in Appendix I for the 
centrosymmetric and non-centrosymmetric cases. The theoretical maximum 
of — 0:89 for the non-centrosymmetric case agrees fairly well with the value 
of — 0-75 as calculated from the actual data used for the anti-correlated 
synthesis. The difference between the two is to be expected since the pro- 
cedure adopted is subject to the same types of errors as are encountered in 
the calculation of a Wilson plot. The synthesis was calculated only for 
the non-centrosymmetric projection and is shown in Fig. 1 (d). It is very 
surprising indeed to find that even this diagram shows the structure in good 
detail. 


3. EXPLANATION OF THE SYNTHESES 


(a) From the Fourier Point of View.—The normalised phase synthesis 
confirms the most essential property of the phase synthesis deduced theo- 
retically in Parts II and IV, namely, that it resembles closely the Fourier synthe- 
sis. The consequence of taking /2'f;? for the coefficients instead of unit 
amplitudes as in e** will be that peaks will not be sharp but will have a finite 
size and shape corresponding to an average atom in the structure. The 
normalisation also makes it possible to compare directly this synthesis with 
the actual Fourier synthesis. It is interesting to see that even the relative 
strengths of carbon, oxygen and chlorine have come out correctly, even though 
the amplitudes fed in corresponc to those of a single average atom in the 
unit cell. 


It will be noticed that the random synthesis is also almost as good as 
the normalised phase synthesis though it contains some spurious peaks as 
compared with the practically uniform background in the latter. That 
the random synthesis would reveal the structure in such good de‘ail could 
not have been anticipated to start with. In fact, it was this rather unexpected 
result thet suggested the next test that was made, namely, to use a set of 
amplitudes completely anti-correlated with the correct ones. The resultant 
synthesis, Fig. 1 (d), is of course worse than the random synthesis, but still 
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the tyrosine molecule is clearly discernible. What is more interesting is 
that the relative strengths of the peaks at the positions of chlorine, carbon, 
etc., are practically unaffected. However, the number of spurious peaks 
is now somewhat larger as compared with the random synthesis and their 
strengths are also larger. 


The real reason why even under these diverse conditions the phase syn- 
thesis is able to reveal the structure can be understood by making a closer 
study of the problem. The explanation follows essentially from the general 
idea of the deconvolution of the Patterson function and of the principle 
of modulation of one structure by another (Parts I and IT). All the above 
syntheses can be interpreted as the modulation of the appropriate modulus 
synthesis with the phase synthesis, e. The latter is practically same as 
the Fourier synthesis and contains first order peaks of strengths 0-90 /;/S, 
at the same positions as in Fourier where Sy = V Zf;? (Part IV). Now, 

j 


because of the fact that the structure amplitudes used in the case of random 
and the anti-correlated syntheses do not correspond to any actual “‘ structure”, 
their modulus synthesis would not contain only positive peaks as in a modulus 
synthesis of any real structure (Part [V). However, at the origin alone both 
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Fic. 2. (a) The fourier synthesis, (b) the normalised phase synthesis and (c) the random 
phase synthesis of L-tyrosine hydrochloride projected down the b-axis. Contours are at intervals 
of 2e/A® at chlorine and 1e/A? elsewhere, starting from 3e/A*. Shaded area corresponds to less 
than 3e/A?, 
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these will have a strength equal to that of the actual modulus synthesis. 
Elsewhere, one would expect a fluctuating background. The only strong 
interactions which arise in the modulation of this synthesis with the phase 
synthesis will therefore correspond to those of an origin peak with the first 
order peaks of the phase synthesis and will thus lead to the structure con- 
tained in the phase synthesis, i.e., at the correct atomic positions. 


Another general feature that might be observed in Figs. | and 2 is that 
the phase synthesis for the non-centrosymmetric case is distinctly better 
than that for the centrosymmetric case. This result is best understood after 
we consider another interesting aspect of the problem, namely, from the 
statistical point of view. 


(b) From the Statistical Point of View-—The terms random and anti- 
correlated have been used in the previous sections in relation to the mean 
of the distribution. These would have conveyed the idea that the distri- 
butions chosen are not correlated or anti-correlated to the distribution of 
the structure amplitudes for the correct structure. However, it should be 
remembered that in a Fourier synthesis it is the actual structure amplitudes 
that is employed so that the correlation about the mean (xk) is not so impor- 
tant as the direct correlation ,c). The terms correlation about the mean 
and the direct correlation are explained in Appendix. 


If Table I of Appendix is examined it will be seen that the direct correla- 
tion between the structure amplitudes of the assumed distribution and the 
actual distribution is as high as 0-79 for the non-centrosymmetric case and 
0-63 for the centrosymmetric case when these are randomly correlated about 
the mean. Consequently, it is not at all surprising that the random synthesis 
reveals much of the details of the actual structure. Even when the two sets 
of structure amplitudes are anti-correlated (or as much anti-correlated as 
possible) the direct correlation coefficient is still as high as 0-60 for the non- 
centrosymmetric case. In this case, even the anti-correlated synthesis does 
teveal the structure quite well and this becomes understandable from the 
magnitude of the direct correlation coefficient. The correlation is poorer 
in the centrosymmetric case (about 0-40) but it is believed it is good enough 
to reveal some details of the structure. 


The large value of the direct correlation arises essentially because the 
quantity | F | is always positive. Consequently, even if there are variations 
about the mean there is always a positive cerrelation between the absolute 
magnitudes in the two sets of the data because both of them must necessarily 
be positive, 
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It will be noticed from Table I that the direct correlation for the non- 
centrosymmetric case is systematically larger than that for the centro- 
symmetric case for all the different types of syntheses. This explains the fact 
already mentioned, namely that the phase synthesis for the non-centrosym- 
metric case is distinctly better than that for the centrosymmetric case, as 
a general feature. 


4. Tests ON A LIGHT-ATOM STRUCTURE 


(a) The Normalised Phase Synthesis.—It will be observed that al! the 
above syntheses have been calculated for the compound L-tyrosine hydro- 
chloride which is not a purely light-atom structure. Hence the theoretical 
considerations given above are not strictly applicable to this case. The 
effect of the presence of heavy atom on the distribution function P (y) has 
already been worked out (Srinivasan, 19605). In the present case, in which 
there are two heavy atoms in the unit cell, the general tendency is to alter 
the distribution curve for both centrosymmetric and non-centrosymmettic 
cases such that there are a large number of reflections having their value 
near the mean. The effect of this in both cases, is to make the value of c 
larger (see Appendix 1) than if the heavy atoms were not present. I[n view 
of this fact, it was felt that the test should be carried out on a pure light-atom 
structure. This was actually done and the normalised phase synthesis for 
both the non-centrosymmetric and centrosymmetric cases were calculated 
for the structure “‘ tyrosine” alone obtained by removing the two chlorine 
atoms from the tyrosine hydrochloride structure. The resultant syntheses 
[Figs. 3 (a) and 3 (b)] will be found to be practicaily as good as Figs. 1! (6) 
and 2(b). It was therefore felt not necessary to calculate the other syntheses. 











Fic. 3. The normalised phase synthesis of L-tyrosine projected down (a) the c-axis and 
(b) the b-axis. Contours are same as in Figs. 1 and 2 respectively. 


(b) Use of the Structure Amplitudes of an Entirely Different Structure. — 
The results of carrying out the different types of syntheses described earlier 
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suggested one other interesting possibility, namely, that of combining the 
structure amplitudes of an actual structure with the phase angles calculated” 
from an entirely different structure. For this purpose, six atoms were 
assumed at the positions marked by crosses in Fig. 4(a) and these atoms 
were assumed to have a suitable scattering factor (actually 1-7 times the 
‘f? value of carbon) falling off with # as fora carbon atom. This was done 
to make the value of ¥ fj? of this structure (which we may call as the correct 


3 
structure A, and denote its structure amplitudes by | F, |) the same as that 
for the assumed structure. Now, the assumed structure (B) was taken to 
be the structure of pure tyrosine. The phase angles og for tyrosine were 
calculated and were used along with the structure amplitudes | F,| to cal- 
culate a synthesis using | F, | exp. iag as Fourier coefficients. The synthesis 
for the non-centrosymmetric projection is shown in Fig. 4(a@). Figure 








FT 


Ng 

















4(b) shows the difference Fourier synthesis (| F,4|— | Fg |) exp. iog. 
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Fic. 4. (a) The synthesis using | F,|exp.aj, as fourier coefficients where A corresponde 
to an arbitrary six-atom structure (marked by crosses) and B to tyrosine. Contours are sams 
as in Fig. 1. (6) The difference synthesis, (| F,|— | F,|) expia,. Contours are at intervals 
of le/A®. Positivs and negative contours are shown by continuous and broken lines respectively. 
The zero contour is omitted. 


5. DISCUSSION 


The synthesis using the structure amplitudes of an arbitrary structure 
along with phases calculated from a different structure represents probably 
the extreme possible case in the present series of tests. The agreement 
between this synthesis [Fig. 4 (a)] and the Fourier [Fig. | (a) with the chlorine 
removed] is very good and it will be seen that all the atoms of tyrosine have 
come out in it. Particularly itis noteworthy that, excepting at one place, 
there are practically no strong peaks at the atomic positions of the arbitrary 
structure A [marked by crosses in Fig. 4(a)], whose structure amplitudes 
were used in the synthesis. The one peak mentioned above is, however, 
believed to be spurious, 
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From this, the following important conclusion is obvious, namely, 
that if a Fourier synthesis is made using the amplitudes of the correct structure 
(4) along with the phases calculated for an assumed wrong structure (B) then 
the resultant synthesis would contain peaks only at the assumed positions 
and not at the correct positions. 


This seems to be in apparent contradiction with the well-known fact 
in structure analysis, namely, that the correct structure is usually arrived 
at by a series of successive Fourier syntheses, assuming initially atoms to 
be present at all the prominent peaks in a trial Fourier synthesis. The 
resolution of this paradox follows from the general idea of modulation 
(Part 1). 


Now, the modulus synthesis, | F, |, contains an origin peak of strength 
0:92Sy and non-origin peaks of strengths 0-38 f, fj/Sw at rai — raj, while 
the phase synthesis, e***, contains peaks of strengths 0:90 f{/Sy at rgi. The 
modulation of the two should obviously lead to peaks at rg; arising out of 
the interactions between the origin peak of the former and the first order 
peaks of e** at rg;. The other interactions between raj — raj and rpi would, 
in general, lead only to a background and in particular there will not be any 
concentration of electron density at r4;. If, on the other hand, the structures 
A and B have some atoms in common, say P of them, then P of the peaks 


in rg; would be of the type raj so that the interactions between raj — raj 
and r,;i will lead to concentration at raj. The strengths of 4 will, of course, 
depend on the number of coincident atoms (P) and will, in general, be weaker 
than those at B. 


It follows from the above that the necessary condition for peaks to occur 
at the correct atomic positions, A, is that at least some of the atoms used for 
calculating the phases must coincide with those in the structure A. 


This again shows that the usual method of successive Fourier syntheses 
would successfully converge to the correct structure only if the assumed 
structure or at least part of it is near the correct one. The method would 
obviously fail if completely wrong positions are used, to start with. Of 
course, there is no reason why the phases fed in from a completely arbitrary 
structure should lead to the correct one since, if it were to be true, it would 
be an absurdly simple solution to the phase problem. 


However, the fact, that the assumed structure is wrong, would easily 
be revealed by the R value. For a completely wrong structure it would 
be 83% for the centrosymmetric and 59% for the non-centrosymmetric case 
respectively (Wilson, 1950). Further, the difference Fourier synthesis might 
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also help in detecting the falsity of the structure. Even this is not so very 
helpful as is usually supposed. Though there is a small concentration of 
negative density in Fig. 4(b) at some of the tyrosine atoms, both positive 
and negative peaks occur at random in the synthesis the maximum values 
corresponding approximately only to a helium atom. Probably the error 
synthesis would have been more useful in detecting the wrong atomic 
positions. 
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APPENDIX I 
Correlation between Structure Amplitudes 


We define the correlation between the structure factois F, and F, of 
two structures | and 2 as 


1 2(F,F,* + F,F,*) _ 2=|F,| | Fy| cos (a, — ay) (1 
2VE(F.P-VElF vzIF PP. vz[RP © 


where the notation is self-explanatory. We would be mainly interested in 
the case when the two phase angles a, and a, are the same for any reflection 
so that the cosine factor in (1) becomes unity giving the following quantity 
which would also be called the correlation between the structure amplitudes 
|F,| and |F,|. Thus 

e( FD = opal Q) 

V2|F,[? V2|F3 |? 

The term coefficient of linear correlation (k) between two statistical variates 
is usually used to denote the correlation about the mean. Applying this 
to our case, we may write k as 


c (F) = 


VE|F, |? V2] F, |? 
We shall hereafter use the term “direct correlation” to denote c(|F}) 
as defined by (2) and the term “ correlation coefficient ” or simply “ correla- 
tion” to denote k as defined by (3). (If it desired to be explicit we may 
use “correlation about the mean” to denote k.) 
It is possible to express ¢ in terms of k by expanding (3) and substituting 
from (2). This gives 
2 m, (| F; |) mm (| Fe |) — ko (| F, |) o (| F2|) (4) 
eu eS = 
Vm (| F, |) V mz (| F, \) 
where o is the standard deviation and m, and mz, are the first and second 
moments of the variates. 
in the discussion to follow we would be mainly interested in the case 
when the two variated have the same probability distributions so that we 
may take m,(|F,|) = (|F.|) = and o (| F, |) = ¢(| F, |) =a, ete. 
Equation (4) then reduces to 


c=k+ m(l — k). (5) 
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When the two variates are one and the same, both k and c should be equal 
to unity which also follows from (5). 


The values of ni, and m, for the normalised structure amplitudes are 
well known in the case of both non-centrosymmetric and centrosymmetric 
distributions. These are m, = (4/z/2), m,=1 for the former and m, = 
(2/7), m, = 1 for the latter. Substituting these values in (5) the values 
of ¢c are obtained to be as in Table I. The values in brackets (where given) 
correspond to the actual values found in the calculation of the correspond- 
ing syntheses described in the text. 


It will be seen from Table I that even for perfect anti-correlation* the 
value of ¢ is quite appreciable and particulaily for the non-centrosymmetric 
case it is quite high. Also the value of c for the non-centrosymmetric case 
is systematically higher than the corresponding value for the centrosymmetric 
case. This essentially arises*because of the fact that the probability distri- 
bution function in the former case has a large number of reflections with 
their values clustered round the mean while in the centrosymmetric case 
they have a comparatively wider dispersion. Thus if we could have an ideal 
case of a distribution which is a delta iunction at the mean, the value ot 
c would always be unity. Table [ also tells us that the worst possible direct 


TABLE I 


Values of the direct correlation c for different values of k 
Actually observed values are given in brackets 








_k 


Soi. 1:0 0-5 0:0 —0-5 —0-80* —0-99* —1.9 
bution 





Non-centro- .. 1:0 0:89 0-79 0-68 43 0-59 0-57 
symmetric (0-74) (0-55) 


Centrosymmetric 1:0 0°82 0-63 O-45 0:34 ee 0:26 
(0.69) 





* These values correspond to the maximum possible values of the anti-correlation coefficient 
for the centrosymmetric and non-centrosymmetric cases respectively. 


* The terms ‘‘correlated’’ and “anti-ccrrelated’’ (atout the mean) are used in general, 
where k is positive and negative respectively, while ‘‘ perfect correlation’’ and “ perfect anti- 
correlation’’ are used to denote the extreme cases of k = +1 and —1 respectively. 
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correlation is for a centrosymmetric structure when the structure amplitudes! 
are perfectly anti-correlated. : 


However, the above discussion on the relationship between c and 
is purely from a theoretical point of view and it remains a separate issu 
whether the case of perfect anti-correlation is at all realisable in practicg! 
for any type of distribution. In fact, it may be stated here that, unless the! 
distribution is symmetrical about the mean, perfect anti-correlation about: 
the mean is not possible. (The details will be presented elsewhere.) Thig 
is also equivalent to the condition that the mean and the median of the diss) 
tribution should be the same. Further, it is possible to calculate the maxi-! 
mum value of anti-correlation attainable in practice for any type of distrie) 
bution. The values for the non-centrosymmetric and the centrosymmetri¢: 
distributions are mentioned in Table I. It will be noticed that this maximum) 
value is larger for the non-centiosymmetric than the centrosymmetric casey 
This is understandable since the distribution for the former case is fairly: 
symmetrical about the mean while in the latter case it is rather asymmetrical, 


Tn conclusion, I wish to express my sincere thanks to Prof. G. N. Rama= 
chandran for valuable suggestions during this investigation. My thanks) 
are also due to Miss V. Valambal and Mr. S. K. Majumdar for helping me’ 
in the calculations. . 
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